Topological Identification of Spin-1/2 Two-Leg Ladder with Four-Spin Ring Exchange 
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states, dimerized Heisenberg models, [231124 125|, and t-J 
model[26]. For these systems, the non-trivial (tt) Berry 
phase on a link reveals a singlet on the link, which is a 
pure quantum object due to two-spin exchange. 

In this paper, we extend the quantized Berry phase 
to be sensitive to the effect of four-spin ring exchange 
interaction, and apply it to a S = 1/2 spin ladder 
with ring exchange interactions. According to the phase 
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diagram[13j, there are two first order transitions from fer- 
romagnetic phase. One is transition to the rung singlet 
phase. The other is that to the dominant collinear spin 
phase, which connect to the dominant vector chirality 
phase through crossover at a self-dual point [17j. These 
two phases are the singlet phases with short-range or- 
der and have dominant correlation of collinear spin and 
vector chirality respectively. The rung singlet phase in- 
cludes the spin ladder with only two-spin exchange inter- 
actions and the dominant vector chirality phase includes 
that with only four-spin exchange interactions. These 
phases have a unique ground state with finite gap, while 
the other phases including ferromagnetic phase do not 
have finite gap under the translational symmetry. More- 
over, through the spin- chirality duality transformation, 
the gap in the dominant vector chirality phase of the 
original Hamiltonian is smoothly connected to that in 
the rung singlet phase of the transformed Hamiltonian 
fiH . Although it is believed that the ground state of the 
rung singlet phase is well approximated by the product 
of local rung singlets, there is no simple picture for the 
dominant vector chirality phase. To clarify it, we discuss 
the adiabatic connection of the Hamiltonian to a simple 
model, i.e., a topologically equivalent model. 



A spin-1/2 two- leg ladder with four-spin ring exchange is studied by quantized Berry phases, 
used as local order parameters. Reflecting local objects, non-trivial (tt) Berry phase is founded 
on a rung for the rung-singlet phase and on a plaquette for the vector-chiral phase. Since the 
quantized Berry phase is topological invariant for gapped systems with the time reversal symmetry, 
topologically identical models can be obtained by the adiabatic modification. The rung-singlet 
phase is adiabatically connected to a decoupled rung-singlet model and the vector-chiral phase is 
connected to a decoupled vector-chiral model. Decoupled models reveals that the local objects are 
a local singlet and a plaquette singlet respectively. 

PACS numbers: 73.43.Nq, 75.10.Jm, 75.40.Cx 



I. INTRODUCTION 

The recent progress in the multiple-spin exchange in- 
teractions is attracting much attention. It has been 
found to be important in several materials such as 
two-leg ladder compound La x Cai 4 _. c Cu2404i 1, H, 
two-dimensional antiferromagnet L^CuO^fj, [5|, mag- 
netism of two-dimensional quantum solids, e.g. solid 3 He 
films [6], and Wigner crystals |7j. Four-spin ring exchange 
plays an essential role in several models to give rise to ex- 
otic phases due to its frustration, e.g., the vector chirality 
phases @, nematic orderings@, and octapolar order [Toj. 

Especially, the two-leg ladder model with the multiple- 
spin exchange interactions has been studied extensively 
El G3, El Q EE EE E3 ■ To clarify its rich phases, not 
only correlation functions corresponding to phases but 
also entanglement concurrence [3 and string order [l3j 
are useful to characterize the phases. As such a novel 
order parameter, which is beyond the Ginzburg-Landau 
symmetry-breaking description, there is an order param- 
eter based on the topological invariants corresponding to 
the topological order [l9j. 

Recently, Berry phases [2(| have been used in order to 
detect the topological order and the quantum order (2li 
[221 . [23( 1 . The Berry phases are quantum quantities based 
on the Berry connection which is defined by the overlap 
between the two states with infinitesimal difference and 
do not have any corresponding classical analogues. Then, 
one can define it even though there is no classical order 
parameter. The advantage of the Berry phase is that 
it quantizes to or tt even in the finite sized systems 
in any dimension when the system has the time reversal 
invariance. It has been successfully applied to several 
quantum systems such as generalized valence bond solid 



II. DEFINITION OF THE BERRY PHASE 
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Let us start with the definition of Berry phase [20( in 
a quantum spin system. For the parameter dependent 
Hamiltonian H(4>), the Berry phase 7 of the ground 
state is defined as 17 = J^ v A((j))d(j) (mod 2vr), where 
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A(4>) is the Abelian Berry connection obtained by the 
single-valued normalized ground state |gs(</>)) of H(<f>) 
as A((f)) = (gs(0) |c?<^ |gs(0)) . This Berry phase is quan- 
tized to or 7r if \gs(cf>)} is a gapped ground state and 
the Hamiltonian H(<fi) is invariant under the anti-unitary 
operation 0, i.e. [7J(</>),0] = 0. It has a remark- 
able property that the Berry phase has topological ro- 
bustness against the small perturbations unless the en- 
ergy gap between the ground state and the first excited 
state closes. We note that the Berry phase is undefined 
if the energy gap vanishes while varying the parameter 
cf>. Then, we limit ourselves to the rung singlet phase 
and the dominant vector chirality (collinear spin) phase, 
which have finite gap. To calculate the Berry phase 
numerically p7| . we use 7 = HmM^oo ~fM , where 7m is 
defined by discretizing the parameter space of <f> into 
M points as j M = - J2m=i argC(0 m ), 4> m = 2nm/M 7 
where C(<fr m ) is defined by C(<j> m ) = (gs(0 m )|gs(0 m+ i)) 
with 4>m+i = 4>\- In general, M can be a very small 
number [281]. 

To study a local structure of the quantum system, 
such as a local singlet, we use a local spin twist on a 
specified link i,j as the parameter </> 23]. Under this lo- 
cal spin twist on a link the term Sj + S^Sj in 

the Hamiltonian is replaced with e^SfSj + e'^S'Sj, 
where s£_ = Sf±iSf. Although previous studies [2ll. I22I 
I23I I25I [26| deal with the spin twist only for the two- 
body terms, we extend the twist to the four-spin ex- 
change interactions. Since the Hamiltonian is written 
by H = J2ij JijSi- Sj + Y^ijki Kijkl(Si-Sj)(S k -Si), the 
local spin twist for a selected link i,j is introduced into 
all the term Si ■ Sj in the Hamiltonian. As described 
below, the extended Berry phase can detect not only the 
local singlet but also the plaquette singlet. 

The quantized Berry phase is considered as a link- 
variable. Then each link is labeled as one of three labels: 
"0-bond" , 'V-bond" , or "undefined" . Especially, we shall 
calculate the leg Berry phase 7/, the rung Berry phase j r , 
and the diagonal Berry phase 7<j. The quantization of the 
Berry phase is guaranteed by the time reversal symmetry 
of the quantum spin system [2l|. 



III. THE MODELS AND THE RESULTS 

A. 5 = 1/2 spin ladder model with four-spin 
exchange interaction 

The S — 1/2 spin ladder model with four-spin ex- 
change interaction is described by the following Hamil- 
tonian 



Hcyc J 



N/2 N/2 



i=l a =l,2 



N/2 
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FIG. 1: (color online) Leg, rung, and diagonal Berry phase 
7i> 7r, 7d of the S = 1/2 periodic TV = 16 ladder model 
with four-spin interactions as a function of 9 ( J = cos#, K = 
sin#) with the schematic pictures of corresponding phases. 
Berry phases are zero, n (shaded region) , or undefined (region 
without data). 



with the ring exchange 

Pi + Pi 1 = Si,i ■ Sia + • Sj+1,2 + SiA ■ Si+i,i 

+ Si,2 ' Si+1,2 + Si ; l ■ + <Sj,2 ' 

+ 4(5^1 • Si t 2)(Si+l t l ■ Si+1.2) 

+ 4(5^1 • Si+\ t i)(Si.2 ■ Si+1.2) 
- 4(S iA ■ S i+h2 )(S it2 ■ S lti+1 ), (2) 

where S^a are the spin-1/2 operators on the site fi, a) 
and N is the total number of sites. The periodic bound- 



ary condition is imposed as S 



N/2+i 



for all of the 



(1) 



models in this paper. We set the parameters as J = cosfl, 
K = sm6. Figure Q] shows Berry phases on local links of 
iV=16 ladder obtained numerically by the exact diago- 
nalization method. Three kinds of Berry phases, 7/ on 
the leg link (i, 1) — (i, 2), j r on the rung link (i, 1) — (i, 2), 
and jd on the rung link (i, 1) — (i + 1, 2), are calculated 
and showed translational symmetry. Nontrivial (jr) Berry 
phases obtained for finite size systems identified two dif- 
ferent phases; One is the rung singlet phase and the other 
is the dominant vector chirality phase including the dom- 
inant collinear spin phase. These three phases have a spin 
gap in thermodynamic limit (ill. [l3| . The crossover be- 
tween the dominant vector chirality phase and the domi- 
nant collinear spin phase is not found. It is neither found 
by the entanglement [IH nor by the Lieb-Schulz-Mattis 
twist operator [l]|. 

Berry phases in the other phases are also obtained in 
Fig. [Tj but becomes undefined due to gap-closing in the 
thermodynamic limit, such as in the ferromagnetic phase 
and around the self-dual point 8 = arctan(l/2) ~ 0.14-7T. 
The staggered dimer and scalar-chirality phases around 
the self-dual point are Z2 symmetry breaking phase with 
two-fold degeneracy in the thermodynamic limit [l8| . In 
a N = 16 system, there is a finite-size gap in the S z = 
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FIG. 2: Minimum gaps for twist <j> on a leg, rung, and di- 
agonal link of the 5 = 1/2 periodic ladder model with four- 
spin interactions as a function of adiabatic parameter a for 
H(a) = (1 — a)H cyc \$ = i,s v +aHns- The several data are 
plotted for N = 8, 12, 16, and 20. 
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FIG. 3: Minimum gaps for twist <f> on a leg, rung, and di- 
agonal link of the 5 = 1/2 periodic ladder model with four- 
spin interactions as a function of adiabatic parameter a for 
H(a) = (1 — a)H cyc |9=o.8vr +aHovc- The several data are 
plotted for N = 8, 12, 16, and 20. 



sector of H cyc for most values of 9 and we can obtain the 
Berry phase. However, results around the self-dual point 
in FigQ] shows the numerical instability, i.e., Berry phase 
becomes undefined. Note that we can define non-Abelian 
Berry phase to avoid the numerical instability with using 
the gap above the two-fold degenerated ground states 
in the thermodynamic limit around the self-dual point. 
The transition point 9 C between the rung singl et p hase 
and Z2 symmetry breaking phase is about 0.l7r[lll [l3|. 
while Fig. Q] shows 9 C < 2ir. To clarify this difference, we 
shall study the Berry phase at 9 — 0. 

Before studying the model at 9 = 0, let us now inter- 
pret the phases by the adiabatic modification in order to 
obtain a decoupled model with the same Berry phases. 
We shall calculate the </> dependence of the energy gap 
to see whether the gap closes or not during the adiabatic 
modification since the Berry phase remains the same if 
the gap does not close during the adiabatic modification. 

a. Rung singlet phase We consider the adiabatic 
modification from i?cyc|e=i.8ir m the rung singlet phase 
to the completely decoupled model which has the Heisen- 
berg type coupling only on the rung bonds: 

N/2 

Hrs = <Sj,i ■ Sj t 2- (3) 

i=l 

These two models are connected by adiabatic parameter 
a as H(a) = (1 — a)H cyc + aH^s- As shown in FiglU 
the minimum gap through the twist <j) does not close. It 
means not only that the gap of ff C yc|e=i.87r adiabatically 
connects to the singlet-triplet gap but also that these 
two models with the same Berry phases are topologically 
identical. This adiabatic connection is consistent with 
the fact that the ground state is well approximated by 
the product of local rung singlets. 

b. Vector chirality phase We also consider the adi- 
abatic modification from i? C yc|e=o.87r in the dominant 



vector-chirality phase to a decoupled model: 



JV/4 

-Hdvc = ^ (£20 x £21,2) • (£21+1,1 x £21+1,2) -(4) 

i=l 



We call its ground state as the "dimerized vector-chiral 
state" since it minimize the local operator (£2^,1 X £21,2)' 

(S , 2i+l,l X S'2i+1,2) = (S^i.l ' S'2i+l,l) {£>2i,2 ' "Sai+l, 2) ~ 

(S 2 i,i • S 2 i+ia) (S 2 i,i ■ S 2l +ia)- This operator is used 
as an order parameter in previous studies and a clas- 
sical spin configuration is corresponds to a 90° spin 
structure [l5|. The grou nd state is the product of pla- 
quette singlet states] 29]. Moreover, through the dual- 
ity transformation [ll|, Hdvc is mapped to the summa- 
tion of S^i,! • 5^,2 + S l 2i+i,i • £21+1,2 — 5 , 2i,i • £24+1,2 — 
£21,2 • £21+1,1- This transformed model is identified as 
the rung singlet phase by the Berry phase. It should 
be emphasized again that the gap with no twist <j) = 
in the dominant vector chirality phase is smoothly con- 
nected that in the rung singlet phase of the transformed 
Hamiltonianpl| . 

Figure [3] shows the minimum energy gap through the 
twist 4>ofH(a) — (l~a)H cyc |e=o.8^ +ai?DVC in the adi- 
abatic deformation. Since the gap does not close during 
the modification in N = 20 system at least, we identify 
that the nontrivial diagonal Berry phase 7<j exhibit the 
decoupled vector-chiral state. We should note that the 
modified Hamiltonian breaks the translational symme- 
try although the original Hamiltonian does not. In other 
words, although the gap does not close by the twist of 
diagonal bond in odd-ih plaquette, it does for the diago- 
nal bond in even-th plaquette. In this sense, the ground 
state of Hcyc|e=o.8u- is locally identical to the vector chiral 
state. 
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phase for > tt/A as rung* phase, which implies the 
limitation of the localized rung-singlet picture and en- 
courages us to use another picture such as the resonating 
valence bond theory [3l|. It is consistent with very strong 
quantum fluctuations introduced on the product singlet 
ground state[32j and with "non-pert urbative" (in J r /Ji) 
behavior obtained numerically |30(. To clarify the phase 
which has no 7r-bond, we need further discussion with 
considering another kind of spin twist, because the twist 
should be corresponds to a local structure of the model. 
It should be emphasized that the rung* phase results 
from the quantum phase transition of the twisted Hamil- 
tonian (<j) = tt) for j r and is not contradicting previous 
studies for untwisted Hamiltonian. 



FIG. 4: (color online) Leg, rung, and diagonal Berry phase 7;, 
j r , 7,; of the S = 1/2 periodic N = 16 ladder model without 
four-spin interactions as a function of 6 (J r = cos#, Ji = sind) 
with the schematic pictures of corresponding phases. Berry 
phases are zero, n (shaded region), or undefined (region with- 
out data). 



B. S — 1/2 spin ladder model without four-spin 
exchange interaction 

Although the rung singlet phase is thought to include 
9 = 0[l3|, Fig. Q] shows that the rung singlet phase does 
not include H cyc \g= {J — l,K = 0). To clarify it, we 
study S = 1/2 spin ladder model without four-spin ex- 
change interaction: 

N/2 N/2 

H\ad = J} y] Sj ; g ■ Sj+l ia + J r Sj,l ■ ^,2(5) 

i=l a=l,2 i=l 

where J\ and J r are parametrized as J\ = smO and 
J r = cos6>, respectively. We consider the antiferromag- 
netic case of0<#<7r/2in this paper to concentrate 
on i?iad|e=7r/4 = #cyc|e=o- Note that Hu d \ e=0 = H RS . 
On the other hand, -ffiadl 0=7172 is the model of decoupled 
two chains, which has gapless excitation in the thermo- 
dynamic limit. 

Figure [4] shows the 6 dependence of the Berry phases 
on local links of H\ a d at N = 16 except for the gap- 
less point 8 = 7r/2. 7d = is trivial because there 
is no diagonal interaction. Although the gap of un- 
twisted Hamiltonian <\> — is smoothly connected [30|. 
the gap of twisted Hamiltonian <f> — ir for 7 r closes at 
9 ~ 0.233-7T < 7r/4 and 7 r changes. We denote the 7,- = 



IV. CONCLUSION 

In conclusion, we have shown that the quantized Berry 
phases is useful to classify the phases of spin chains with 
four-spin interaction. In the dominant vector-chirality 
phase which comes from the ring exchange interaction, 
the Berry phase on each diagonal link is used as a pla- 
quette variable and becomes a non-trivial value (tt), while 
the Berry phase has been used as a link variable in pre- 
vious study. The Berry phase is also useful to clarify the 
phase boundary from the finite-size systems since it is 
quantized even in the finite size systems. The property 
of the phase is revealed through the adiabatic deforma- 
tion into a decoupled model: the rung singlet phase (and 
the vector chirality phase) is corresponds to a product 
of rung singlets (plaquette singlets). These two phases 
are connected through the duality transformation. The 
other phases of this model will be detected by the Berry 
phase with another kind of twist or the non-Abelian 
Berry phase for two-fold degenerated ground states. 
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